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Evolutionary games on the lattice: 
death-birth updating process 

Stephen Evilsizor and Nicolas Lanchier* 

Abstract This paper is concerned with the death-birth updating process. This model is 
an example of a spatial game in which players located on the d-dimensional integer lattice 
are characterized by one of two possible strategies and update their strategy at rate one 
by mimicking one of their neighbors chosen at random with a probability proportional to the 
neighbor’s payoff. To understand the role of space in the form of local interactions, the process 
is compared with its nonspatial deterministic counterpart for well-mixing populations, which 
is described by the replicator equation. To begin with, we prove that, provided the range 
of the interactions is sufficiently large, both strategies coexist on the lattice for a parameter 
region where the replicator equation also exhibits coexistence. Then, we identify parameter 
regions in which there is a dominant strategy that always wins on the lattice whereas the 
replicator equation displays either coexistence or bistability. Finally, we show that, for the 
one-dimensional nearest neighbor system and in the parameter region corresponding to the 
prisoner’s dilemma game, cooperators can win on the lattice whereas defectors always win 
in well-mixing populations, thus showing that space favors cooperation. In particular, several 
parameter regions where the spatial and nonspatial models disagree are identified. 


1. Introduction 

This paper is concerned with a closely related version of the death-birth updating process in evo¬ 
lutionary game theory introduced in m- This model is an example of a spatial evolutionary game 
designed based on the framework of interacting particle systems. The concept of evolutionary game 
theory is an extension of traditional game theory that has been proposed by m to describe the 
dynamics of populations in which fitness is frequency dependent: individuals are viewed as players 
who are characterized by their strategy and receive a certain payoff through their interactions with 
other individuals. The payoff is then interpreted as fitness so that players with a larger payoff have 
a higher reproductive success. One of the most popular models of evolutionary game is the repli¬ 
cator equation m which assumes that the population is well-mixing. In contrast, the death-birth 
updating process includes a spatial structure in the form of local interactions: players are located 
on the set of vertices of a connected graph and interact with a finite set of neighbors, meaning that 
they update their strategy based on these neighbors. See m chapter 9] for a general definition 
and a brief review of such models. The model in m assumes that the updates are neutral with 
high probability and based on the payoff of the neighbors with small probability, which we refer 
respectively as voter and game steps. In contrast, the model considered in this paper only accounts 
for game steps, so the duality techniques [6] developed for voter model perturbations are no longer 
available tools to study the process. The main objective is to study the limiting behavior of the spa¬ 
tial stochastic process and confront our results with the limiting behavior of the replicator equation 
in order to understand the effects of the inclusion of space. 
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Model description. The process studied in this paper, which we again refer to as the death-birth 
updating process following the terminology in El, is a spin system on the d-dimensional integer 
lattice where each vertex is occupied by a player characterized by one of two possible strategies, 
say strategy 1 and strategy 2. The state at time t is a function 

^ {1, 2} where ^tix) = strategy at vertex x at time t. 

The dynamics of this process or any other spatial game is defined in a couple of steps: we first fix 
a payoff matrix, which allows us to turn every spatial configuration of strategies into a so-called 
payoff landscape, which can then be used to define the transition rates at each vertex. Since we 
focus on games with two strategies, the payoff matrix is a 2 x 2 matrix A = (a^) whose coefficients 
are positive real numbers interpreted as 

ttij = payoff of a type i player interacting with a type j player. 

In nonspatial evolutionary games, players interact equally with any other player in the population, 
making their payoff a function of the global frequency of representatives of each strategy. In contrast, 
spatial games assume that the payoff of a player depends exclusively on the strategy of a finite set of 
neighbors, which is the key to designing more realistic models with local interactions. Throughout 
this paper, the interaction neighborhood of vertex x is the set 

Nx -■= {y y ^ X and maxj=i^2,...,d \xj -yj\ < M] 

where the constant M is called the range of the interactions. Letting Nj{x,^) be the number of 
neighbors of the player at vertex x following strategy j, every spatial configuration ^ is then turned 
into a payoff landscape by attributing the payoff 

4>{x,0 ■=Ylj aij Nj{x,^)l{C{x) =i} for i = l,2, (1) 

to the player at vertex x. In words, each type i player receives aij from each of her neighbors 
following strategy j. The last step to define the dynamics of the process is to follow m and 
interpret the payoff as fitness. The basic idea here is to write the rate at which a player changes her 
strategy as a function of her payoff and the payoff of her neighbors in such a way that players with 
a larger payoff are more likely to spread their strategy. There are multiple options. For instance, 
the updating rules considered in [9l [13] are as follows. 

• Best-response dynamics [9|. Players update their strategy at rate one in order to maximize 
their payoff, which depends on the strategy of their neighbors. 

• Payoff affecting birth and death rates [T3] . In this process, when a player has a positive payoff, 
at rate this payoff, one of her neighbors chosen at random adopts her strategy, whereas when 
her payoff is negative, at rate minus this payoff, she adopts the strategy of one of her neighbors 
chosen at random. This updating rule is inspired from [3]. 

The dynamics of the death-birth updating process is built using a similar approach: we assume that 
players update their strategy at rate one by mimicking a random neighbor, with each neighbor being 
chosen with a probability proportional to her payoff. More precisely, letting ^ be the configuration 
of the system, the player at x switches her strategy i ^ j at rate 


Pi-^j{x,0 ■ 


for {i,j} = {l,2}. 


( 2 ) 
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This model, or to be more specific, a closely related version of this model, has been introduced 
and studied heuristically in m while HE] give rigorous results. The process considered in these 
works can be seen as the weak selection approximation of the model described by ([2|). Players again 
update their strategy at rate one but, at the time of the update, 

• with probability 1 — e, the player mimics the strategy of a neighbor chosen uniformly at 
random, just like in the voter model mm, while 

• with probability e, the player mimics a neighbor chosen at random according to probabilities 
that are proportional to the neighbors’ payoff, as described by ([2|). 

This model is studied in HElttZ] when e is small, in which case duality techniques for voter model 
perturbations are available. In contrast, we study the process when e = 1, in which case duality 
cannot be used, which leads to more qualitative and less quantitative results. 

The replicator equation. Before studying the spatial game, it is worth taking a quick look at 
its nonspatial deterministic analog to later identify disagreements between both models and thus 
understand the effect of the inclusion of space in the form of local interactions. The nonspatial model 
is obtained by assuming that the population of players is well-mixing, which results in a system 
of ordinary differential equations for the frequency of each strategy. In the case of the death-birth 
process, this is a time-change of the replicator equation: 

u[ = Ui U2 {4>i{ui,U2) - 4>2{ui,U2)) (3) 

where Uj is the frequency of players following strategy j and 

4>iiui,U2) = aiiui + ai 2 U 2 and (/> 2 (mi,U 2 ) = a 2 i'Ui-|- 022^2 

are the common payoffs of all type 1 and all type 2 players, respectively. This can be viewed as the 
nonspatial analog of the payoff landscape ©• As pointed out in [T2|, the limiting behavior can be 
conveniently described by introducing the parameters 

:= ail ~ 0,21 and 02 := 022 — ai 2 

and calling strategy i selfish whenever a* > 0 and altruistic whenever Oj < 0. Then, following the 
usual terminology by calling a strategy an evolutionary stable strategy if it cannot be invaded by 
any alternative strategy starting at an infinitesimally small frequency, some basic algebra shows 
that the behavior of the replicator equation ([3|) is as follows: 

• when ai 02 < 0, the selfish strategy always outcompetes the altruistic strategy, showing that 
the selfish strategy is the only evolutionary stable strategy, 

• when 01,02 > 0, there is an unstable interior fixed point so the system is bistable, showing 
that the two (selfish) strategies are evolutionary stable, 

• when 01,02 < 0, there is a globally stable interior fixed point so both strategies coexist and 
none of the two (altruistic) strategies is evolutionary stable. 

In summary, the analysis of the replicator equation shows that, when the population is well-mixing, 
a strategy is an evolutionary stable strategy if it is selfish but not if it is altruistic. 

Main results for the spatial game. In order to compare the spatial game with its nonspatial 
analog, we assume that the process starts from a spatially homogeneous distribution, i.e., a product 
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measure in which the density of each of the two strategies is constant across the lattice. Since the 
two configurations in which all players follow the same strategy are absorbing states, we also assume, 
to avoid trivialities, that the density of each strategy is positive. For the spatial game, 

• strategy j wins when limt^oo P = j) = 1 , 

• strategy j survives when liminfi_^oo P = j) > 0 , 

a strategy is said to go extinct when it does not survive and both strategies are said to coexist when 
they both survive. Note that the two probabilities above do not depend on the choice of vertex x 
because both the initial distribution and the evolution rules are translation invariant. From now 
on, we assume without loss of generality that 021 > 012 > 0 and study the limiting behavior of the 
process as the other two payoffs vary. 

To begin with, we look at the parameter region where both strategies are altruistic, in which case 
coexistence occurs when the population is well-mixing. Numerical simulations suggest that, except 
in the one-dimensional nearest neighbor case, coexistence is again possible for the spatial game 
though the coexistence region is reduced. The smaller the spatial dimension and/or the range of 
the interactions, the smaller the coexistence region. Our first two theorems show that coexistence 
is indeed possible and that the coexistence region for the spatial game is indeed reduced. More 
precisely. Theorem [ 1 ] shows that, regardless of the spatial dimension, both strategies coexist when 
they are sufficiently altruistic and the range of the interactions is sufficiently large. 

Theorem 1 - There is a > 0 such that coexistence occurs when 

max (ail, 022) < cl and M is sufficiently large. 

To prove that the coexistence region is reduced, and more generally identify parameter regions 
in which one strategy wins, we first observe that, when on = 012 and 022 = 021) the process is 
significantly simplified because the payoff of the players only depends on their strategy but not on 
the strategy of their neighbors. Indeed, in this case ([T|) reduces to 

Nj{x,^) l{C(x) =i} = an (( 2 M -F 1 )'^ - 1 ) l{C(x) = i} 
for i = 1 , 2 , therefore the transition rates become 

ccjj Nj{x,C) 






aiiNi{x,i) + Nj{x,C) 


for {i,j} = { 1 , 2 }. It follows that, under our general assumption 012 < 021, the set of type 2 players 
dominates stochastically a certain biased voter model BIS], thus showing that, in this very special 
case, strategy 2 wins. Elaborating on this idea but using coupling arguments to compare the death- 
birth updating process with spin systems which are more complicated than the biased voter model, 
we can prove much more, as shown in the next theorem. 

Theorem 2 - Recall that 021 > ai2 and let N := cardNx- Then, 

(a) strategy 1 wins when min (012 — 022, uii — 021) > (-N — l)(a2i — 012), 

(b) strategy 2 wins when (M, d) / ( 1 , 1 ) and 

{N'^ - N - 1 ) max (ail - 021, ai2 - «22, an - CL22) < a2i - ai2. 
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Figured] shows the parameter regions in both parts of the theorem. Note that the parameter region 
in the hrst part of the theorem is nonempty if and only if 

ai2 — {N — l)(a2i — 012) >0 if and only if 012 > (1 — 1 /-^) 021- 

The hgure shows this region for N = 2 , i.e., in the presence of one-dimensional nearest neighbor 
interactions. In contrast, the parameter region in the second part of the theorem is always nonempty 
and, more interestingly, it always overlaps the region where both strategies are altruistic as well as 
the region where both strategies are selhsh. This shows that the inclusion of a spatial structure in 
the form of local interactions indeed reduces the coexistence region, as mentioned above. This also 
shows that, in a subset of the parameter region where the replicator equation is bistable, there is 
instead a strong type for the spatial game that wins even starting at low density. The next theorem 
goes a little bit further in this direction by showing that, no matter how selfish a strategy is, the 
other strategy always wins if it is selfish enough. 

Theorem 3 - For all a > 0 there is A < 00 such that strategy 1 wins when 

max (021,022) < a and an > A. 

Note that, contrary to Theorem [ 2 l this theorem does not assume that 012 < 021, therefore it also 
extends the parameter region where strategy 2 wins found before. 

The results collected so far indicate interesting discrepancies between the death-birth updating 
process and the replicator equation, showing the importance of local interactions. The most inter¬ 
esting aspect suggested by spatial simulations is the existence of a subset of the parameter region 
corresponding to the prisoner’s dilemma game in which cooperators win on the lattice whereas they 
always lose when the population is well-mixing. The prisoner’s dilemma game is characterized by 
the following ordering and terminology of the four payoffs; 

ai2 = sucker’s payoff < 022 = punishment 

< ail = reward < 021 = temptation. 

Figure [ 2 ] shows the corresponding triangular region in solid lines. Players with strategy 1 are called 
cooperators while players with strategy 2 are called defectors. Because the reward is not as good as 
the temptation, and the punishment is not as bad as the sucker’s payoff, cooperators are altruistic 
and defectors selhsh, therefore defectors indeed win when the population is well-mixing. In contrast, 
the heuristic arguments in m suggest that there is a subset of the prisoner’s dilemma triangle in 
which cooperators are favored over defectors on regular graphs. This has been proved in [ 3 ] for 
hnite, connected, simple graphs, and in [ 6 ] for integer lattices with d > 2 . Their results, however, 
hold in the weak selection case but not for the process ([2]). We now study the interactions among 
cooperators and defectors in one dimension, the main difficulty being the lack of attractiveness of 
the process. To state our last result, we introduce the following quantities that will be interpreted 
later as drift of a certain interface: 


D3 := 


L »4 := 


Oil + «12 


«21 + O22 


an -|- ai2 + 0,21 -\- (122 2 oii -|- 021 -|- 022 


«11 + 012 


021 + 022 


On -f 012 -|- 2a22 
Then, we have the following theorem. 


2 an -|- a 2 i -|- 022 
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Figure 1 . Phase diagram of the spatial game along with a summary of the theorems in the on — 022 plane. In the 
picture, the points p- and p+ are the two points introduced in the proof of Lemma 


Theorem 4 - Assume M = d = 1 . Then, strategy 1 wins when 

(022 < £121 and D3 + D4 > 0) or (022 > 021 and D4 > 0). 


Note that the parameter region given by the theorem overlaps but is not restricted to the prisoner’s 
dilemma triangle. To see that the theorem implies the existence of a subset of the triangle in which 
cooperators win, observe that, when an = 021 > 022 = cli2, 


1 , £112 + £121 „ £112 + £121 ^ 1 ,1 „ 1 n 

D‘i + Da = - - 2 X - >-H- 2 x-= 0 . 

2 3ai2 + 021 £1 i 2 + 3a2i 2 2 2 


In particular, the first parameter region given by the theorem in which cooperators win indeed 
overlap the prisoner’s dilemma triangle. The theorem also implies that strategy 2 wins in the 
parameter regions obtained by exchanging the role of the two strategies. Noticing in addition the 
symmetry in the expression of D4, we deduce from the second part of the theorem that 


strategy 1 wins when 022 > £121 and D4 > 0 
strategy 2 wins when an > 042 and D4 < 0 


showing that, when min(an,a22) > max(ai2, 021), the condition is sharp. Figure [ 2 ] gives a picture 
of the curves derived from the theorem when a2i/ai2 = 2. 
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“22 “11 = “12 “11 = “21 



Figure 2. Phase diagram of the one-dimensional nearest neighbor spatial game when 021/012 = 2. The curves are 
obtained from the conditions in Theorem [4] and their symmetric expressions. The triangle in solid lines represents 
the parameter region corresponding to the prisoner’s dilemma game. 


2. Coexistence of altruistic strategies 

This section is devoted to the proof of our coexistence result Theorem [TJ For simplicity, we focus on 
the two-dimensional case but our approach easily extends to any spatial dimensions. Specifically, 
we will prove that both strategies coexist when M is large and 

max (an, 022) < 5 “^ 2 “^^ (c-)^ min (ai2,021) = 2 “^"^ {c+)~^ min (ai2, 021) ( 4 ) 

where the two key constants c_ and c+ are defined as 

c_ := 2 “^^ min (012/021,021/012) and c+ := 5 ^ 2 ’^(c_)“®. ( 5 ) 

Let s := ln( 2 ) and := [—rM, rM)^ for all r > 0 , and fix 

A,Bc i^i/2 with card(^) = card(S) = 2 “^M^. 

The proofs of Lemmas [6HZ1 below hold for such general sets though they will be applied ultimately 
to more specific space-time boxes. One key to the proof is to observe that 

maxj=i^2 \xj — Vjl < M for all {x,y) G Ax B. 
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For all Del? finite and i = 1 , 2 , we let 

Ct{D) := card{x G D : ^t{x) = i} 
denote the number of type i players in the set D at time t. 

Keeping the players in a box - To begin with, we prove in the next lemma that, the number 
of players of either type in a given spatial region does not decrease too fast. The idea is to simply 
find a bound for the number of updates using standard large deviation estimates for the binomial 
random variable. This lemma will be used repeatedly later. 

Lemma 5 - For all D d'l? finite, n G N and i = 1 , 2 , 

P{Q{D) < K for some t G ( 0 ,ns) | C,l{D) > K) < exp(- 2 -("+ 3 ) K). 

Proof. To begin with, we let 

Ui{D) := card{x G D : ^o{x) = i and Ti{x) < ns} 

denote the total number of players in the set D who are initially of type i and update their 
strategy at least once by time ns. Recalling that the random variables Ti{x) are independent and 
exponentially distributed with rate one, our choice of s implies that 

Ui{D) = Binomial {(q{D), 1 — = Binomial {Q{D), 1 — 2 “""). 

In particular, using the large deviation estimate 

P (Binomial {K,p) < K {p — z)) < exp{—Kz'^/ 2 p) for all 2: G ( 0 ,p) ( 6 ) 

with p = 2 “"' and 2 = we get 

P {Q{D) < K for some time t G ( 0 ,ns) | QiP) P K) 

< P {ui{D) > (1 - 2 -(”+i)) K I Ci{D) = K) 

= P (Binomial (RT, 1 - 2 -”) > (1 - 2 -(’"+i)) K) 

= P (Binomial (RT, 2 “”) < 2 -(’"+i) K) < exp(- 2 -(’"+ 3 ) RT). 

This completes the proof of the lemma. □ 

Moving the players around - We now prove that if the region A has a large number of type 1 
players then, regardless of the conhguration around this region, we can “move” a positive fraction 
of these players to the nearby region B in s units of time. The constant c_ defined in ([ 5 ]) will play 
the role of the fraction of players we can move. 

Lemma 6 - Assume (j3|) and let c_ as in ([ 5 ]) and a > 0 . Then, 

P(Cs(B) < c-{aM) I Co(^) > < exp(—(aM)^/^) for all M large. 
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Proof. Applying Lemma [5] with n = 1 and K = aM, we get 

P{CHB) < C-{aM) I CoH^) > aM and CqH-B) > aM) 

^ PiCtiP) ^ 2~‘^{aM) for some t € (0,s) | Co(^) > aM and Cq{B) > aM) 

< exp(—2“^aM) < ex.p{—{aM)^P) 

for all M large. It remains to prove that 

P{Cl{B) < C-{aM) I C^{B) <aM< Co'(^)) < exp(-(aM)V2) ( 7 ) 

for all M large. To lighten the notation, we let 

P*(E) := P (E I Co(P) <aM < CqH^)) for any event E 
and introduce the two events 

A := < 2~‘^ aM for some time t G (0,s)} 

B := {(f{B) < for some time t G (0, s)}. ^ ^ 

Applying Lemma [5] with K = aM then K = 2~^ M^, we get 

P*(A) < exp(-2-%M) and P*(B) < exp(-2-^ A/^). (9) 

In addition, observing that the conditional payoff of each type 1 player in the set A given that the 
number of type 2 players in the set B is large satisfies 

(l){x I Ct{x) = 1, B'^) > ai 2 2~^ for all t G (0, s) 

we deduce that, on the event (A U B)'^, each time a player in the set B updates her strategy, she 
remains/becomes of type 1 with probability at least 

Pi > iai 2 2-5 M2)(2-2 aM)((ai 2 2-5 m2)(2-2 aM) + (aai + a22)(2M + 

> (ai2 2-5 M2)(2-2 aM)(a2i 25 M^)-^ > 2-^2 ( 012 / 021 ) oM-^ (10) 

> 25 c_(oM-i) 

for all M large. In particular, letting 

Xu ■= card |x G P : Ti(x) < s| 

V ; j 

Ai := card {x € B : Ti{x) < s and ^s{x) = 1} 
observing from (fTOl) that 2-5 M'^ pi > c-{aM), and using ([ 6 |), we deduce that 
P(Ai < c-{aM) I (AUB)^) 

< P{Xu< 2-4 m2) + P (Ai < c_(aM) | (A U B)‘= and Xu > 2-4 m2) 

< P (Binomial (M2/4, 1/2) < 2-4 M^) + P (Binomial (2-4 M'^,pi) < c-{aM)) (12) 

< P (Binomial (M2/4, 1/2) < 2-4 M^) + P (Binomial (2-4 M'^,pi) < 2-5 M'^pi) 

< exp(—2-5 m 2)-|-exp(—2-’^ M2pi) < (1/2) exp(—(aM)4/2). 
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Figure 3 . Picture related to the proof of Lemma 0 


Using that Xi < Cg{B) and combining ([9]) and (fT^ . we conclude that 

P*{CKB) < c_{aM)) < P*(A U B) + P (Ci(P) < c_(aM) | (A U B)'^) 

< P* (A) + P* (B) + P (Ai < c_ (aM) | (A U B)'^) 

< exp(— 2 “^aM) + exp(— 2 “'^ M^) + ( 1 / 2 ) exp(—(aM)^/^) < exp(—(aM)^/^) 
for all M large, which gives ([7]). This completes the proof. □ 

Creating a pile of players - The next lemma improves the previous one by showing that if the 
region A has a large number of type 1 players then the same amount of type 1 players can be 
created in the nearby region B. The idea is to first prove that, as long as the number of type 1 
players nearby is small, the number of such players can be increased by a factor c+. Once this 
threshold is reached, one can find a small box with a large number of type 1 players and apply the 
previous lemma repeatedly to move a fraction of these players to the target set B. 

Lemma 7 - Assume that 0 holds. Then, 

P (C 6 s(P) < M I Co(^) > M) < exp{—M^P) for all M large. 

Proof. To keep track of the amount of type 1 and type 2 players in A 5 / 2 , which is the key to 
controlling the payoff of the type 1 players in the set A, we set 

K := {Cti.P 5 / 2 ) P c^M for some time t G (0, s)} 
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where c+ has been dehned in ([5]). The proof is divided into two steps. 

Step 1 - First, we prove that the event K occurs with probability close to one. The idea is to show 
that the complement of this event confers a large payoff to the type 1 players in A, which results in 
a large production of such players with high probability, thus leading to a contradiction. To make 
this argument precise, we observe that, with probability one, 

(^(x + Ki) > for all x G ^ and Ct{^ 5 / 2 ) < c+M 

on the event and for all t G ( 0 , s). 

4>{x I ^t{x) = 1, K'^) > ai 2 for all x € A 

4>{x I = 2, K'^) < 021 c+M -|- 022 (2M -|- 1)^ < 2^ 022 M^ for all x G K^/ 2 - 

In particular, given (A' U K)^ where A' is the first event ([ 8 |) for o = 1, each time a player in B 

updates her strategy, she remains/becomes of type 1 with probability at least 

qi > (oi 2 M 2 )( 2 - 2 m)((oi 2 M 2 )( 2 - 2 m) + ( 23 o 22 M 2 )( 2 M + 1 ) 2 )-i 

> (oi2M2)(2-2M)((24a22M2)(2M + l)2)-i (13) 

> 2 “^ (012/022) M“h 

Defining Xu and Xi as in ([TT]) . observing that, by (jl|) and (fT3ll . 

2“^ M^ > 2"^^ ( 012 / 022 ) M > c+M 
and using ([ 6 ]), we deduce that, for all M sufficiently large, 

P (Ai < c+M I (A' U K)‘=) 

< P{Xu< 2“^ m2) + P{Xi< c+M I (A' U K)'= and Xu > 2“^ M^) 

< P (Binomial (M^/d, 1/2) < 2“^ M^) -|- P (Binomial (2“'^ M^, qi) < c+M) (14) 

< P (Binomial (M^/d, 1/2) < 2“^ M^) -|- P (Binomial (2“^ M^, qi) < 2~^ M^ q^i) 

< exp(—2 “® m2 )-|-exp(—2 “'^ m2 gi) < (1/4) exp(—M^/2(_ 

Using again Ai < Cs{B) and combining Q with o = 1 and (fTT)l . we get 

P(K'= I CoH^) >M) = P{CHB) < c+M and | CoH^) > M) 

< P (A' I Co'(^) > M) + P (C/(P) < c+M I (A' U K)'^) 

< exp(—2“^ M) -|- (1/4) exp(—M^/2^ 

< ( 1 / 2 ) exp(—M^/ 2 ) 

for all dispersal range M sufficiently large. 

Step 2 - Now, given K, there is a box with at least 5 2 c+M type 1 players. One can move a 
fraction of these players to the target set B in at most hve steps, applying Lemma [H To begin 
with, we observe that some basic geometry implies that, given K, there exist 

Po, • • • ,-85 C A 5/2 and toG{0,s) 

such that the following three conditions hold: 
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(a) We have Cto(-Bo) > 5“^c+M. 

(b) For A; = 0,..., 5, we have card(i?fc) = 2“^ with = B. 

(c) For A; = 0,...,4, we have maxj=i ^2 \xj — yj\ < M for all (x,y) £ B^ x i?fc+i- 

We refer to Figure [3] for an illustration of the worst case scenario where all the type 1 players are 
located in one of the corners of K^i 2 . Under these conditions, we can bring type 1 players to our 
target set in at most five steps applying repeatedly Lemma El Indeed, 

P iChP) < for all t G (0, 6s) | K) 

< PiCl+5sm < 2^M I dm > 5-2 C+M = 27(c_)-5m) 

< l-nfc=o,i,2,3,4^(Ci(i^fc+i)>2^(c-)"-"M|Co'W>27(c-)'=-''M) 

< 1 - n.=0,l,2,3,4 (1 - exp(-2V2 (,_)(.-5)/2 ^1/2)) (16) 

< 1 - (1-exp(-27/2MV2))5 

< 5 X exp(—2^/2 < (1/4) exp(—M^/2^ 

for all dispersal range M sufficiently large. 

Conclusion - Combining (|15p - (ll6l) . we deduce that 

P < 2^M for all t G (0, 6s) | CoH^) > M) 

< pm \ Co(^) >M)+P iCdP) < 2^M for all t G (0, 6s) | K) 

< (1/2) exp(-MV2) ^ ( 1 / 4 ) exp(-MV2) ^ ( 3 / 4 ) exp(-MV2) 

which, applying Lemma El with n = 6 and K = 2^M, implies that 

PiCi{B)<M\m)>M) 

< P{d{B) < 2^M for all t G (0,6s) | Co(^) > M) 

+ P {d{B) < 2-’^2^M I CtiB) > 2^M for some t G (0,6s)) 

< (3/4) exp(—M^/2) _|_ exp(—2-2 M) < exp{—M^P) 

for all M large. This completes the proof. □ 

Block construction - To deduce coexistence, we use Lemma[7]in combination with some obvious 
symmetry and a block construction. The idea of the block construction is to define a coupling 
between the process under consideration properly rescaled in space and time and supercritical 
percolation. More precisely, let be the directed graph with vertex set 

H := {{z, n) X ■. zi + Z 2 + n \s even} 

and in which there is an oriented edge 

(z, n) —>■ ( 2 :', n') if and only if {z' = z P ei ov z' = z P 62 ) and n' = n + 1 

where ej is the jth unit vector. Then, we consider the 14 dependent oriented site percolation process 
with density equal to 1 — e on this directed graph, i.e., we assume that 

P {{zi,ni) is closed for i = 1, 2,... , m) = e™' 
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whenever jzi — Zj\\/ \ni — nj\ > 14 for i ^ j. Then, we set 

:= (M/2) 2 -h [-M/4, M/Af for all z G and T := 6 s = 6 ln(2) 

and declare site (z, n) G H to be good whenever 

Cnri^z) = card {x G Bz '■ inT{x) = i} > M for i = 1, 2. 

Finally, for all n G N, we define 

:= {z : {z,n) is wet} and := {z : {z,n) is good}. 

The next lemma shows that, for all e > 0, one can find a sufficiently large dispersal range such that 
the set of good sites dominates stochastically the set of wet sites. In view of the dehnition of a good 
site, this will imply coexistence of both strategies. 

Lemma 8 - Assume (j4]) and fix e > 0. Then, for all M large, 

P {z G W^) < P {z G Xn) for all {z, n) G H whenever Wq G1 Xq. 

Proof. First, we define the collection of events 

'Bi{z,n) := {(l^rp(^Bz) > M} for all {z,n) G H and f = 1,2. 


Then, since for j = 1,2, 

Bz, Bz±ej c ((M/ 2 ) 2 ; ± (M/4) Cj) + and card(.B 2 ) = caxd{Bz±ej) = 2“^ M^ 

we can apply Lemma [7] to get 

P (Bi( 2 ; ± Cj, n -|- 1) for j = 1, 2 | 61 ( 2 , n)) 

= > M for i = 1,2 I CnTiB^) > M) 

> 1 - 4 X P > M I > M) > 1 - 4 X exp(-MV2) 

> l-e/2 

for all M large. By symmetry, the same holds for strategy 2, therefore 

P{z± ej G Xn+i for j = 1,2 I 2 ; G Xn) 

= P (Bi( 2 ; ± Cj, n -|- 1) n B 2 ( 2 ; ± e^, n -|- 1) for j = 1,2 | Bi( 2 ;, n) n B 2 ( 2 ;, n)) 

> — 1-|-P (Bi( 2 ; ± Cj, n-|-1) for J = 1,2 I Bi( 2 :, n) n B 2 ( 2 ;, n)) (17) 

+ P (B 2 ( 2 ; ± Cj, n + 1) for j = 1,2 | Bi( 2 ;, n) n B 2 ( 2 ;, n)) 

> - 1 + 2(1 - e/2) = 1 - e. 

Since in addition all the estimates in Lemma [7] hold regardless of the conhguration outside the 
spatial region /2 j we deduce from (fT7)l the existence of a collection of events 

0 ( 2 ;, n) for all {z,n) G H 


that satisfy the following three properties: 
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(a) The event G{z,n) is measurable with respect to the graphical representation in 

R{z,n) := {{M/2)z,nT) + {Ky^ x [0,T]). 

(b) For all M large, we have P (G(z, n)) > 1 — e. 

(c) We have the inclusion G{z, n) fl {z G X^} C {z ± ej G Xn+i for j = 1, 2}. 

Observing also that 

R{z, n) n R{z', n') = 0 when |z — z'| V |n — n'| > 2 x 7 = 14, 

we deduce from [ 8 l Theorem 4.3] the existence of a coupling between the long range death-birth 
process and the oriented site percolation process such that 

P {W^ C Xn) = 1 whenever Wq C Xq. 

The lemma directly follows from the existence of this coupling. □ 

To conclude, we fix e > 0 small enough to make the percolation process supercritical. Observ¬ 
ing that, starting from a product measure with a positive density of both strategies, the number of 
good sites at level zero is almost surely infinite, we deduce that 

liminfo^oo -P (0 G W2n \ Wo = Xq) > 0. 

This, together with Lemma [HI implies that, for all M large, 

liminft^oo P {Ct{x) = l)Pfe(x) = 2) 

= liminft^oo P (6(0) = 1) P (6(0) = 2) 

> liminft^oo P (6(0) = i | 0 G kF 2 p/rj) P (0 G kF 2 p/rj | Wo = Xq) > 0 

for all X G . This completes the proof of Theorem [1] 

3. Coupling with modified voter models 

This section is devoted to the proof of Theorem [2j The common ingredient to prove both parts of 
the theorem is to couple the process with the modified voter models Ct whose transitions 

at vertex x are given by the following expressions 

i^j at rate Ci^j{x,C) ■= i'i-- e) fj{x,C) + e l{fi{x,C) = 0} 
j ^ i at rate Cj^i{x, C) := (1 - e) fi{x, C) + e l{fi(x, C) 7 ^ 0 }. 

for {i,j} = { 1 , 2 } and where 

fj{x,C) = cardjy G : C*(a^) = j}/cardX^, = {l/N)Nj{x,C) 

denotes the fraction of neighbors of vertex x in state j. In words, the transition rates indicate that 
particles are updated at rate one and that, at the time of an update, 
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• with probability 1 — e > 0, the new type is chosen uniformly at random from the interaction 
neighborhood just like in the voter model, 

• with probability e > 0, the new type is i unless all the neighbors are of type j. 

The results of |12l section 3] show using duality that type i particles win for this process. In 
particular, to prove that strategy 1 wins, it suffices to prove that the set of type 1 players in the 
death-birth process dominates stochastically its counterpart in Ql, which follows from 

^ and ^(x) = implies that 

Pi^ 2 {x,i) < ci-, 2 {x,C^) and p 2 -,iix,C) > C 2 ^i(x,C^) 


according to Theorem III.1.5 in M- Since in addition the transition rates of the modihed voter 
models are monotone with respect to the number of neighbors of each type, in order to show that 
strategy 1 wins, it suffices to prove that the simplified implication 


A^i(x,.^) = A^i(x,C^) and ^(x) = C^(^) implies that 

Pi-, 2 ix,C) < ci-, 2 {x,C^) and p 2 -,i{x,C) > C 2 -^i(x,C^) 

holds for some e > 0. By symmetry, strategy 2 wins if the implication 

^ 2 {x ,0 = N 2 {x,C'^) and ^(x) = C^(x) implies that 

Pi-, 2 ix,C) > ci-, 2 {x,C‘^) and p 2 ^i(x ,0 < C 2 -^i(x,C^) 


(18) 


(19) 


holds for some e > 0. Using (|18p . we now prove Theorem [2ja. 


Lemma 9 - Recall that 021 > 012 . Then, strategy 1 wins when 

min (ai 2 - 022 , oii - 021 ) > (A^ - l)(o 2 i - 012 ). 

Proof. In view of the discussion above, it suffices to prove that (fTSl) holds. First, we observe that, 
when the fraction of type 1 neighbors of vertex x is equal to either zero or one, the transition rates 
are the same for both processes so it remains to prove (|18p under the assumption 

A''iA '"2 / 0 where Nj := Nj{x,^) = Nj{x,C^) for j = 1 , 2 . ( 20 ) 

The transition rate at vertex x depends on the payoff of its neighbors, and the main idea is to 
express the transition rates by distinguishing between the part of the payoff coming from x and the 
part of the payoff coming from the other neighbors’ neighbors. In order to make this distinction, 
we introduce the following four weighting factors: 

Wij ■■= T.y^x (lU(y) = 0 Y.zr^y,z^x 1{?(^) = J}) for bJ = 1>2. 

That is, Wij is the number of type j neighbors (excluding vertex x) of a type i neighbor of vertex x 
counted with order of multiplicity. Note that, for z = 1,2, we have 

+ Ej=i,2 Wij = Ni + J2y^x (iU(y) = i} J2z^y Hz ^ 3 :}) 

= N, + {N-l) Ey^^ my) = i} 

= Ni + (N-1)N, = NN,. 


(21) 
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In addition, for i ^ j, the transition rates can be expressed as 




_(-/Vj’ + Wji) aji + Wjj ttjj _ 

(A ^2 '^ii) “ 1 “ '^ij ^ij ^ji “t“ '^jj ^jj 


Using (l 2 T]) - (l 2 ^ . we now prove that (fTSl) holds in the nontrivial case (l 20 t) . 
Transition 1 —)• 2 - Using (f 2 T]l and 022 < ai 2 < 02 I) we get 


{N2 + in2l) 021 + W22 022 < {N2 + tC21 + W22) O21 = NN2 021 
(IVi + iRii) On + W12 012 = (A^i + ten + wu) on + wu (012 - on) 

= NNi ail + W12 (012 - Oil). 


This, together with ()22p for i = 1 and j = 2, implies that 

/ ^ _ NN2 021 _ 

~ A^A^i Oil + iri 2 (012 — Oil) + A^A ^2 021 

_ NN2 021 _ NN2 021 

NNi (oil - 021) + tci2 (012 - Oil) + 021 A^2 o2i + Pi 

where, since on — 021 > {N — l)(a2i — 012), 


Pi := NNi (on - 021) + W12 (012 - on) 

= (IViVi - u;i2)(oii - 021) + ini2 (012 - 021) 

> {N - l){NNi - wi2){a2i - 012) + W12 (012 - 021) 

= N {{N - 1 ) Ni - Wi 2 ){a 2 i - 012) = N wii (021 - 012) > 0 . 

Note that the strict inequality holds because (l2n]l - (f2T]l imply that 


NNi — W 12 = (A^i + wii + u;i 2 ) — W 12 = Ni + wu > Ni > 0. 


(22) 


(23) 


(24) 


In view of (j 23 l) - (l 24 p . there exists ei > 0 small such that 

P 1 ^ 2 (X ,0 < IViV 2 02 l(IVVl+Pl)-l = A^ 2 (^+(Pi/A' 02 i ))-1 

< (1 - e) f2{x, 0 = (1 - e) f2ix, C^) = ci^ 2 (a;, C^) 

whenever e < ei and (f20P holds. 

Transition 2 — > 1 - Since on > 021 > 012 and 022 — 012 < {N — l)(ai2 — 021), using the previous 
estimates and obvious symmetry, we show that 

P2^i(a;,0 >-^^1 012 (N^^oi 2 + P2)~^ where p2 < N W22 (au - a2i) < 0 . 


In particular, there exists 62 > 0 small such that 

P2-,iix,C) > Ni{N + {p2/Nai2))~^ 

> (1 - e)/i(o ;,0 + e = (1 - e)/i(x,C^) + e = C 2 ^i(x,C^) 

whenever e < 62 and (f 20 P holds. 

In conclusion, the implication (1181) holds for e smaller than min ( 61 , 62 ) > 0, which shows that 
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strategy 1 wins under the assumptions of the lemma. □ 

Repeating the proof of Lemma 0 step by step but exchanging the role of both strategies only 
shows that strategy 2 wins under the strong assumption an < ai2 < 021 < 022- In fact, this 
sufficient condition for strategy 2 to win can be easily improved to 

max (an, 012) < min (021, a22) 

by using a coupling between the death-birth process and a biased voter model with a selective 
advantage for type 2 particles to show that the former dominates the latter. To prove that strategy 2 
wins in the larger region stated in Theorem [ 2 jb, we couple the death-birth process with the second 
modihed voter model but using techniques different from the ones used to show the first part 
of the theorem. To explain the assumptions of the theorem, we note that our approach requires 
the interaction neighborhood to have a certain connectivity property which does not hold in the 
one-dimensional nearest neighbor case. First, we introduce the payoff functions 

4 >i{z) := an {z/N) + au (1 - z/N) = (an - ai2)iz/N) + an 

(j)2iz) := 022 (z/N) + 021 (1 - z/N) = (022 - a2i){z/N) + 021 

for all 2: = 0,1,... , Ai, and let a := max (an, 012,021,022) and 

M_|_ := max^ (/i{z) > max^^iv </'i(^) =: 
m_ := min^ (l)2{z) < min^^jv 4‘2{z) ='■ rn+. 

The following lemma gives a sufficient condition on these minimum and maximum payoffs for 
strategy 2 to win. This condition is made more explicit in the subsequent lemma. 

Lemma 10 - Strategy 2 wins whenever {M,d) 7 ^ (1,1) and 

(N — 1) > {N — 2 ) M_|_ -|- M_ and {N — 1) M_ < {N — 2 ) m_ -|- m_|_. ( 25 ) 

Proof. Following as in the proof of Lemma O it suffices to show that (jl 9 p holds. This is again 
trivial when the fraction of type 1 neighbors of vertex x is equal to either zero or one so we focus 
from now on on the nontrivial case where 

N1N2 7^ 0 where Nj := Nj{x,^) = Nj{x,(‘^) for j = 1 , 2 , ( 26 ) 

indicating that x has two neighbors y* and z* with different strategies. Except in the one-dimensional 
nearest neighbor case M = d = 1, one can find two vertices y* and 2* such that 

y*, 2* 7^ X and y* G Nx n Ny* and 2* G H and y* G 

and we may assume that y* and 2* are neighbors of each other: 

C{y*) = 1 a-iid >^(2:*) = 2 and y*,z* G and y* G Nz*. 

The rest of the proof is divided into two steps depending on the transition. 

Transition 1 —)■ 2 ~ Assume that ^{x) = C‘^{x) = 1. Then, 

4>iy I C{y) = 1) < M+ and 4 >{y \ ^(y) = 2) > m+ for all y G N^. 
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In addition, we have (p{y*,C) < M- therefore 

Pi^2{x,i) > 


m+ N 2 


{Ni — 1) M+ + M_ + N 2 
m+ N 2 


( 27 ) 


(m+ - M+) N2 + {N - 1) M+ + M_ 
Now, for all 2 ; G [0, A^], we define the functions 


■ (m+-M+)z + (Af-l)M+ + M_ 

Since (A^ - 1) m+ > (A^ - 2) M+ + M_, 

m+ 


1 , / N (1 “ f) -2^ 

and hi{z) := --he. 


S-iCl) = 


> 


m+ + (iV - 2) M+ + M_ N' 

We then distinguish two cases. 

Case 1 - When M+ > m+, it follows from (l28]i that 

gi{z) > gi{l) z > z/N for all zG[1,A^ —1] 

therefore gi{z) > hi{z) for all e > 0 small by continuity. 

Case 2 - When M+ < m+, the function gi is concave down. Moreover, 


(28) 


for all e > 0 small according to (f28ll . while 

giiN) = 


Nm. 


> 1 = hi{N). 


Nnij^ — M+ + M_ 

This again implies that gi dominates hi for all e > 0 small (see Figured for a picture). 
Recalling ([271) . we deduce that, in both cases and when (l26l) holds, 

Pi^2{x,0 > 51 (^ 2 ) > hi{N2) = ci^2ix,C‘^) 

which proves the first inequality in (I19p . 

Transition 2 —)• 1 - Assume that ^{x) = C‘^{x) = 2. Then, 

(t>{y I i{y) = 1) < M_ and ct){y \ i{y) = 2) > m_ for all y G N^. 

In addition, we have (/>(z*,^) > m+ therefore 

M Ni 


P2^l{x,i) < 


M_ Ni + {N 2 — 1) m_ + m+ 

M_Ni 

(M_ — m_) A^i + (A" — I) m_ + m+ 


( 29 ) 
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Now, for all z G [0, A^], we define the functions 

M_z 

g2{z) := 


(M_ — m_) z + {N — 1) m_ -|- m+ 
Since {N — 1) M_ < {N — 2) m_ -|- < (N" — 1) 


and h 2 {z) 


52(1) 


Af_ 1 

-|- (A^ — 2) rri- + N 

iN-l)M. 1 

(A^ - 1)M_ -hm+ A^' 


As before, we distinguish two cases. 

Case 1 " When M_ > m_, it follows from (1301) that 


N 


(30) 


92 {z) < g 2 {'^) z < z/N for all zG[1 ,A^ —1] 

therefore g 2 {z) < h 2 {z) for all e > 0 small by continuity. 

Case 2 - When M_ < m_, the function g 2 is concave up. Moreover, 

g 2 {N — 1) > (1 — 1/A^)(1 — e) = h 2 {N — 1) for e > 0 small 

52(0) = /i2(0) = 0 

where the first inequality follows from (I30p . This again implies that /i 2 dominates g 2 for all e > 0 
small, and we refer to the right-hand side of Figure [5] for a picture. 

Recalling (12^ . we deduce that, in both cases and when (l26l) holds, 

P2^i{x,i) < g2{Ni) < h2{Ni) = C2-^i(x,C^) 

which proves the second inequality in (jlOp . 

This completes the proof. □ 


To complete the proof of Theorem [2l the last step is to re-express the condition in the previous 
lemma using the payoff coefficients. 

Lemma 11 - Assume that 012 < 021 - Then (1251) holds whenever 

(A^^ - A^ - 1) max (an - 021,012 - 022 , an - 022 ) < 021 - 012 - (31) 

Proof. We distinguish four cases depending on the sign of an — ai 2 and a 22 — 0 , 21 - 
Case 1 - When an < ai 2 and a 22 > 021 , we have 

A/+ = A/_ = ai 2 and m_ = m+ = a 2 i 

therefore (j25l) holds if and only if 012 < a 2 i, which is true by assumption. 

Case 2 - When an > ai 2 and a 22 > a 2 i, we have 

M+ = an, A/_ = an-I-( 1 /A^)(ai 2 - an), m_=m+ = a 2 i. 
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Figure 4. Picture related to the proof of Lemma fTOl 


Using some basic algebra, we deduce that (l 2 ^ holds if and only if 

— N — l)(aii — 021) < 021 “ 012 and {N — l)(aii — 021) < 021 “ 012 
therefore (f 25 ]) holds if and only if {N'^ — N — l)(aii — 021) < 021 — 012- 

Case 3 - Assume that on < 012 and 022 < 021. This case can be deduced from the previous one 
by symmetry exchanging the role of the two strategies, and we find that 

(| 25 ]l holds if and only if — N — l)(oi2 — 022) < 021 — 012. 

Case 4 ~ When on > 012 and 022 < 021, it is easier to prove the result graphically and we refer 
to the phase diagram of Figured] for an illustration of some of the arguments of the proof. In this 
case, the minimum and maximum payoffs are given by 

Af+ = On and M_ = on + (1/A')(ai2 — on) 
m_ = 022 and m+= 022 + (l/A^)(o2i - 022) 
so the two inequalities in ( 1251 ) are respectively equivalent to 

{N - 1)2 022 > (N^-N-l) an + 012 - {N - 1) 021 
{N - 1)2 an < {N^-N- 1 ) 022 + 021 - {N - 1) 012. 

Since in addition on > 012 and 022 < 021, this specifies two triangles with two common sides, one 
vertical side and one horizontal side that intersect at point 

p:= (012,021) in the On - 022 plane. 

For the hrst inequality in (I 32 p . the third side of the triangle is the segment line going through 
point and with slope s+ where 

p+ := (021 + - N - l)“i (021 - 012), 021) 

s+ := (Af 2 - AT - 1 ){N - l )-2 = 1 + (iV - 2 ){N - l )-2 > 1 . 


( 33 ) 
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Using some obvious symmetry, one finds that the third side of the triangle specified by the second 
inequality in (1321) is characterized by the point and slope 

p-■.= {ai 2 ,ai 2 - {N'^ - N-l)~^ {a 2 i - ai 2 )) and s_ = l/s+< 1. (34) 

Since the segment line connecting p_ and has slope one, the triangle (p,p_,p+) is contained 
in the intersection of the two triangles specified by (j33|) - (|34p . In particular, whenever the payoff 
coefficients are in this triangle, which is equivalent to 

ail > ai 2 and 022 < 021 and - l)(aii - 022 ) < a 2 i - 012 , 

the two inequalities in (1321) hold. 

Since all four cases hold simultaneously when (|3T]) holds, the proof is complete. □ 

Theorem [2jb directly follows from Lemmas [10] and [TT] 


4. Coupling with a pure growth process 

This section is devoted to the proof of Theorem [3| which states that, the other payoffs being fixed, 
strategy 1 wins whenever the payoff coefficient an is sufficiently large. The intuition behind this 
result, which is also the first step of the proof, is to observe that, in the limit as an goes to 
infinity, type 1 players with at least one neighbor of their own type never change their strategy. 
In particular, the set of type 1 players dominates a Richardson model [18], i.e., a contact process 
with no death, which obviously implies that strategy 1 wins. Using a block construction and the 
fact that the transition rates are continuous functions of the payoffs, we deduce that the process 
reaches an equilibrium with a density of type 1 close to one when an is finite but large. The rest 
of the proof consists in showing that we can indeed convert the remaining type 2 players, which 
directly follows from percolation results already established in mm- 

To turn our sketch into a rigorous proof, we let Ct be the d-dimensional Richardson model with 
parameter p. In this spin system, each vertex of the d-dimensional integer lattice is either empty 
or occupied by a particle. Each particle produces a new particle which is then sent to a neighbor 
chosen uniformly at random at rate p. This results in either an empty site becoming empty or two 
particles coalescing in case the target site is already occupied. In addition, occupied sites remain 
occupied forever. More formally, using the usual notation 0 for empty and 1 for occupied, the 
transition rates of the process at vertex x are given by 

co^i{x,C) = and ci^o(a;, C) = 0. 

In order to compare the process properly rescaled in space and time with oriented site percolation, 
we also introduce the space-time regions 

Bk ■■= [-K, kY and Bk{z) := Kz + Bk for all z € Z'^. 

Then, we have the following lemma, where the processes under consideration have been identified 
to the set of vertices in state I to lighten the expressions. 
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Lemma 12 - For all a, e > 0 there exist A, K,c < oo such that 

P {B2K't- for some t G {cK, 2cK) \ Bk C ^0) < e 
whenever max (021,022) < a and an > A. 

Proof. To begin with, we note that 

4'{x,C) P 0:11 > CLii when x G ^ and/i(x, 7 ^ 0 

< (2M + max ( 021 , 022 ) when 

from which it follows that, as an —>■ 00, 

Pi^2(a;,0 ^ 0 when /i(x,O/0 

P2^iix,^) ^ I when /i(x,^/2(a;,0 / 0 for some y G iVa,. 


(35) 


Looking now at the Richardson model, we observe that, because there is no death and because each 
particle newly created must be in the neighborhood of its parent, the set of occupied sites satisfies 
the following connectivity property: 


X G Ct implies that fi{x,^t) ^ 0 for all times t 


(36) 


provided this holds at time zero. Combining (I35]) ~ (p6]l and using that the set Bk is connected, we 
deduce that the death-birth process can be coupled with the Richardson model with parameter one 
in such a way that, for all hxed K,c> 0, 


P {B 2 K n Ct B 2 K n Ct for some t G {cK, 2cK) | Co = Co = Bk) 0 


(37) 


as an —> 00 . In addition, it directly follows from the shape theorem |18] for the Richardson model 
that there exists a positive constant c > 0 such that 


P {B 2 K iC Ct for some t G {cK, 2cK) \ Bk C Co) 

= P {B2K iC CcK I Bk C Co) < -P {B2K iC C,cK \ Co = {0}) < e /2 


(38) 


for all K large. Now, hx RT, c > 0 such that (|38p holds. Since the transition rates of the death-birth 
updating process are continuous with respect to the payoff coefficients and since the space-time 
region in the event in (j37l) is finite, there is A < 00 such that 


p {B2K n Ct B2K n Ct for some t G {cK, 2cK) | Co = Co = Bk) < e /2 
for all an > A. Combining pSil - pUl) . we conclude that 


( 39 ) 


P {B 2 K iC Ct for some t G {cK, 2cK) \ Bk C Co) 

< P {B2K iC Ct for some t G {cK, 2cK) \ Bk C Co) 

-h P {B 2 K n Ct B 2 K n Ct for some t G {cK, 2cK) \ Co = Co = Bk) 

< e/2 -|- e/2 = e 


for all an > A. This completes the proof. □ 
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From the lemma, we deduce that, starting from a product measure with a positive density of 
type 1 players, the density of type 1 at equilibrium is close to one when an is large. To prove this, 
we consider as previously the directed graph ^ with vertex set 

H := {{z, n) € X : zi + Z 2 + ■■■ + Zd + n is even} 
and in which there is an edge (z, n) —>■ {z^, n') if and only if 

z' = z ^ ej for some j = 1,2,... ,d and n' = n + 1. 

Then, calling {z, n) € H an occupied site when 

Bk{z) C for all t G cnK + {0,cK) 

it follows from Lemma da that, for all e > 0, one can choose an large enough so that the set of 
occupied sites dominates stochastically the set of wet sites in the percolation process where sites 
are closed with probability e. Since the probability e can be made arbitrarily small, the density of 
type 1 players at equilibrium can be made arbitrarily close to one. 

The last step is to turn the remaining type 2 players into type 1 players. To do this, the basic idea 
is to rely on the lack of percolation of the dry (not wet) sites for a certain oriented site percolation 
process where sites are closed with a sufficiently small probability e. The fact that the set of dry 
sites does not percolate for small positive e is proved in [3 section 3] for the percolation process 
described above. This result, however, is not sufficient to conclude because the lack of percolation 
of the dry sites for this percolation process does not imply extinction of strategy 2. To solve the 
problem, we consider oriented site percolation on a directed graph that has the same vertex 
set as before but additional arrows, namely 

{z,n) —>■ {z',n') if and only if 

z' = z ± ej for some j = 1,2,... ,d and n' = n + 1 
or z' = z ± 2ej for some j = 1, 2 ,... , d and n' = n. 

process on has the following two key properties: 

As for the process on M’, the dry sites do not percolate if sites are closed with a small enough 
probability e > 0. This is proved in [121 section 3] following the ideas in [7]. 

Recalling that the death-birth process and the percolation process on ^ are coupled in such 
a way that the set of occupied sites dominates the set of wet sites, if 

it{x) = 2 for some ix,t) G Bk{z) x {cnK,c{n + 1)K) 

then site {z,n) can be reached by a directed path of dry sites embedded in This second 
property is also established in [121 section 3]. Even though the proof applies to another model, 
it easily extends to the death-birth process because it only relies on the fact that a type 2 
player can only appear in the neighborhood of a type 2 player. 

To deduce extinction of strategy 2, we first fix e > 0 small such that the set of dry sites does not 
percolate for the percolation process on the directed graph Then, we take an large enough 
so that the set of occupied sites dominates the set of wet sites in the percolation process on the 
smaller directed graph Finally, it follows from the second property above that, because the dry 
sites do not percolate, the type 2 players do not survive. 


The 

1 . 

2 . 
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5. The prisoner’s dilemma in one dimension 


This section is devoted to the proof of Theorem [J] which focuses on the one-dimensional death-birth 
process with nearest neighbor interactions. First, we explain how and D 4 in the statement of 
the theorem are obtained. To do so, we let 

Pi{ni,n 2 ) '■= the rate at which a type i player with 

one type 1 neighbor that has ni type 1 neighbors and 
one type 2 neighbor that has n 2 type 2 neighbors 
update her strategy 

for i = 1,2 and ni, n 2 = 0,1, 2, and note that 


. , (2 - ^2) 021 + 02 O22 

Pl{ni,n2) = -—-r-—^- 

ni On + [ 2 -ni) 012 + [2 - n2) 021 02 022 

ni On-F (2 - ni)ai2 

P 2 {ni,n 2 ) = -—-^-—-r-^- 

ni Oil -I- (2 - ni) 012 + (2 - 02 ) 021 -|- 02 022 
Now, for the process starting with only 2s to the left of the origin, we let 


Xt := inf {x G Z : Ctix) = 2} — 1 and Kt := inf {x > 0 : ^t(x -|- Xt) = 1} 


be the position of the rightmost type 1 player with only type 1 players to her left and the distance 
between this type 1 player and the closest type 1 player to her right. Then, letting 

■= lim/i-^o h~^ E {Xt+h - Xt\it and Kt = j) 

we have the following almost sure estimates 


• o • X X D2{Ct) = D2 = 2 - pi{ 2 , 0 ) 

••• i i i o o • X = D3 = p2{l,l) - pi{ 2 ,l) 

••• i i i o o o X = D4 = P2{1,2) - pi{ 2 ,l) 

for all j > 3 and where x means type 1 or type 2. In particular, D 3 and D 4 are possible drifts of 
the interface at Xt depending on its distance to the next type 1 player. Before studying the process 
starting from general initial configurations, we look at the process starting from configurations that 
have a finite interval of type 1 players, only type 2 players to the left of this interval and infinitely 
many players of each type to the right. In picture, this looks like 


2222111112 
oooo**««*o X X X 

For the process starting from this configuration, we let 

:= inf {x G Z : ^i(x) = 1} and Mt := inf {x > 0 : ^t(x -|- Z^) = 2} 

be the position of the leftmost type 1 player and the distance between this type 1 player and the 
closest type 2 player to her right. Then, we have the following result. 
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Lemma 13 - Assume that 022 < 021 and > 0. Then, 

P {Mt > 3 for alH > 0 and Mt ^ oo \ Mq > 3) > c for some c > 0. 

Proof. To begin with, we let Zt := where a > 0 and 

^j{a) := hm/i_^o h~^ E {Zf+h - Zt\^t and Mt > 3 and Kt = j) 

for all j > 1. To study these quantities, we let 

Z^ := Z^ + Mt — 1 and Kt := inf {x > 0 : ^t(x -|- Xt) = 1} 

be the right boundary of the type 1 interval starting at Z^ and the distance between this type 1 
player and the closest type 1 player to her right. Then, we have 

hm/i_^o P {Z^+h = + 11 6 and Mt > 3 and Kt = 3) = p 2 (l, 1) 
hm/i_^o h~^ P = ■2't^ - 11 6 and Mt > 3 and Kt = 3) = pi(2,1) 


with probability one, while for the left boundary Z^ , 

lim,,^o h-^ P {Zi+h = + 11 6 and Mt > 3) = pi( 2 , 1 ) 

limfe_^o h~^ P {Z~^^ = Z^ -l\^t and Mt > 3) = p 2 {l, 2 ). 

Since Mt = Z^ — Z^ + 1, it follows from (l4n]l - (|4T]) that 

^>3(0) ;= hmfe_^o h~^ E {Zt+h - | 6 and Mt > 3 and Kt = 3 ) 

= Zt (e““ - 1) lim/i_^o h~^ P {Mt+h - Mt = l\ Mt> 3 and Kt = 3 ) 

Zt (e“ - 1 ) lim/j^o P {Mt+h - Mt = -l\Mt> 3 and Kt = 3 ) 
= (e-“ - l)(p2(l, 1 ) + P2(l, 2)) Z* + 2 (e“ - 1 ) pi(2, 1 ) 

therefore, taking the derivative at a = 0, we get 

‘^3(0) = “(P 2 (l, 1 ) +P 2 ( 1 , 2 )) Zt + 2 pi( 2 , 1 ) Zt = —(II3 + Di) Zt. 

Using the same approach, we prove in general that 

“ ~{PjAA + D4) Zt < “(II3 + Pa) Zt < 0 for all j > 1 

since D3 < < D2 when 022 < 021- In particular, 


^jib) < ^j(O) = 0 for some 6 > 0 fixed from now on 

showing that, as long as Mt > 3, the process Zt is a supermartingale with respect to the natural 
filtration of the death-birth process. To conclude, we now apply the optional stopping theorem to 
this super martingale using the stopping times 

r 3 := inf {t : Mt < 3} and := inf {t : Mt > n} for all n > 2. 


26 


S. Evilsizor and N. L angrier 


Since Tn := min(r 3 , Tn) is finite, whenever Mq > 3, 
e-4^ > E{Zo)>E{ZtJ 

> E {Zt^ I Tn = rs) P {Tn = T 3 ) + E {Zt^ I Tn = Tn) P {Tn = Tn) (42) 

> e-3^ (1 - P {Tn = Tn)) + P {Tn = Tn). 

Since the event {Tn = Tn} is nonincreasing with respect to n for the inclusion, we also deduce from 
the monotone convergence theorem that 

P {Mt > 3 for all t > 0 and Mt —>■ oo) 

^ ^ (43) 

> P{Tn = Tn for all n> 3) = hm„_,.oo P {Tn = r„). 

Combining ()12]) - (I43]1 . we deduce that 

P {Mt > 3 for all t > 0 and Mt —>■ oo | Mq > 3) 

> lim„^oo P {Tn = Tn) > limn-s-oo (e“^^ - e“‘^^)(e“^'' - e~^’^)~^ 

> (e-3^ - e-'^^) = 1 - e"'’ 

therefore the lemma holds for c := 1 — >0. □ 

To deal with the process starting from product measures, we note that every realization induces a 
partition of the space-time universe into type 1 and type 2 connected components. More precisely, 
assuming that there is initially a type 2 at the origin, we define the type 2 connected component 
starting at the origin as 

Co := {{x,t) € Z X M+ : there is a path (0,0) —>■2 {x,t) going forward} 
where (0,0) —>-2 {x, t) means that there are times and vertices 

0 = < ^2 < • • • < tn+l = t and 0 = Xi,X 2 , ■ ■ ■ ,Xn = X 

such that the following condition is satisfied: 


{^s{xj) = 2 for all tj < s < tj+i) holds for j = 1, 2,... , n. 


Then, we have the following lemma. 

Lemma 14 - Assume that 022 < 021 and D 3 + D 4 , > 0. Then, 

T := inf {t > 0 : Cq fl (Z x {t, 00)) = 0} < 00 with probability one. 

Proof. We proceed by contradiction, showing that when A := {T = 00} occurs, its complement 
occurs with probability one. To begin with, note that 

• on the event A, the type 2 connected component Cq is unbounded and 

• since there are infinitely many type 1 players on both sides of the origin at time zero, this 
property remains true at all times. 
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Mf It n M+ 



Figure 5. Picture related to the proof of Lemma [HI 


From these two observations, we deduce that 

0 < card {x € Z : (x, t) € Cq} < oo for all t £ (0, oo) 


which, in turn, implies that the left boundary It and right boundary rt of the type 2 connected 
component satisfy the following properties at all times: 

—oo < It := inf {x G Z : (x, t) £ Cq} 

< sup{x G Z : (x,t) G Cq} =: rt < oo. 


We also observe that 

Mr := inf {x > 0 : — x) = 2} > 1 

:= inf {x > 0 : ^i(rt -|- x) = 2} > 1. 

Figure [5] gives an illustration of these processes. Now, the evolution rules of the death-birth updating 
process clearly imply that there exist ci, C 2 >0 such that 


P (min > 3 | min (M^ ,M^) > 1) > ci 

P {T <t + l \rt-lt <'i) > C 2 


(45) 


while it follows from Lemma [13] that 


P (inf {t > s : min (M^ , ) < 3} 

> inf {t > s : r* — < 3} | min (Mp, M+) > 3) > > 0. 


( 46 ) 
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Combining (j44ll - (f46]l . we conclude that 

A = {T = 00} occurs implies that 

implies that 
implies that 
implies that 

This completes the proof. □ 

Lemma 15 - Assume that 022 > 021 and D 4 > 0. Then, 

T := inf {t > 0 : Cq H (Z x {t, 00)) = 0} < 00 with probability one. 

Proof. First, we note that the conclusion of Lemma fT3] holds as well under the assumptions of the 
present lemma, since whenever 022 > 021 we have 

D4 < ZI3 < D2 and Dj^4 + D4 > 2D4 for all j > 1. 

In particular, the lemma follows by repeating the proof of Lemma fT^ □ 

Under the assumptions of Theorem [4] and conditional on the origin being initially 
a type 2 player, it follows from Lemmas [T 3 HIS] that the type 2 connected starting at 
bounded with probability one. Since in addition the number of players is countable, 

P (strategy 2 wins) = P (Co(3^) = 2 and = 00 for some a: € Z) 

— Z]x:5oP)=2 ^ ~ ~ 

where Tx denotes the time at which the connected component starting at vertex x 
defined in the statement of Lemmas 11411151 This proves Theorem [H 
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